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By expressing the time-independent Schrodinger equation in one dimension as a system of two 
first-order differential equations, the transfer matrix for a rectangular potential barrier is obtained 
making use of the matrix exponential. It is shown that the transfer matrix allows one to find the 
bound states and the quasinormal modes. A similar treatment for the one-dimensional propagation 
of electromagnetic waves in a homogeneous medium is also presented. 
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Expresando la ecuacion de Schrodinger independiente del tiempo en una dimension como un sistema 
de dos ecuaciones diferenciales de primer orden, se obtiene la matriz de transferencia para una 
barrera de potencial rectangular haciendo uso de la exponencial de matrices. Se muestra que la 
matriz de transferencia permite hallar los estados ligados y los modos cuasinor males. Se presenta 
tambien un tratamiento similar para la propagacion unidimensional de ondas electromagneticas en 
un medio homogeneo. 
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1. Introduction 



A standard problem in elementary quantum mechanics is that of finding the reflection and trans- 
mission amplitudes for the scattering produced by a potential barrier, or well, in one dimension 
(see, e.g., Refs. 1-4). The reflection and transmission amplitudes are conveniently arranged in the 
transfer matrix, which relates the wave function at both sides of the potential barrier, in such a 
way that the effect of two or more potential barriers is readily obtained by means of the product 
of the corresponding transfer matrices (see, e.g., Ref. 5 and the references cited therein). A similar 
result applies for the one-dimensional propagation of electromagnetic waves in layered media (see, 
e.g., Ref. 6). In fact, the transfer matrices can be defined in all cases where there is an output that 
depends linearly on an input; some important examples, apart from the two already mentioned, 
are the electric circuits and optical systems. In the cases considered here, the transfer matrices are 
2x2 complex matrices but, depending on the equations involved (more specifically, the number of 
variables and the differential order), the size of the transfer matrices do vary. 

The aim of this paper is to show that the transfer matrix for a rectangular potential barrier 
(and, therefore, for a piecewise constant potential) can be easily obtained integrating the time- 
independent Schrodinger equation in one dimension by means of the matrix exponential. The 
time-independent Schrodinger equation in one dimension, being a second-order ordinary differential 
equation, is equivalent to a system of two coupled first-order differential equations and, only in 
the case of a (piecewise) constant potential, this system can be easily integrated using the matrix 
exponential. We also show that making use of the transfer matrix one can find the bound states and 
the quasinormal modes. The transfer matrix for the one-dimensional propagation of electromagnetic 
waves in a medium with a piecewise constant refractive index is obtained in a similar manner, 
without employing the Fresnel coefficients. 

In Sec. 2 an elementary discussion about the transfer matrices for the one-dimensional Schrodinger 
equation is given (see also Ref. 5 and the references cited therein). In Sec. 3 the transfer matrix 
for a rectangular barrier is obtained making use of the matrix exponential; the bound states and 
quasinormal modes are then found starting from the transfer matrix. In Sec. 4 a similar derivation 
for the case of the one-dimensional propagation of electromagnetic waves in layered media is given. 
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2. Transfer matrices 



The solutions of the time-independent Schrodinger equation 

fi 2 d 2 ib 

with a given short-range potential V(x), which vanishes outside the interval a < x < b, can be 



expressed in the form 



^ ie ifc(x-a) + ^-ikix-a)^ for x < Qj 
Bie ik(x~b) + B2e -ik(x-b)^ for x > ^ 

n(x), for a < x < 6, 



(2) 



where = \/2mE/h, A\, A%, B\, B2 are constants and u(x) is a function that depends on the 
explicit form of the potential V{x). By imposing the usual conditions of continuity of ip{x) and its 
derivative at x = a and x = b, a linear relation of the form 

''.ill I 15 1 I 

(3) 
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can be obtained, where M is some 2x2 complex matrix (the transfer matrix), which depends on 
V{x) and the value of k. 

Assuming that V{x) is real, Eq. (fTJ) implies that the probability current density 

2im V dx dx 

where * denotes complex conjugation, satisfies the continuity equation, dj/dx = 0, that is, j(x) 
const.; then, making use of Eq. ([2]), one finds that, for k real 



l^il 2 - \A 2 \ 2 = iSil 2 - \B 2 \ 2 . 



(4) 



Using the fact that 

/ A x W 1 \l A 1 

a 2 ) \ o -1 ) y A 2 

where the ' denotes the Hermitian adjoint, and Eq. ([3]) one finds that Eq. (jl]) is equivalent to 
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V 
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(5) 



The complex 2x2 matrices satisfying Eq. ([5|) form a group with the usual matrix multiplication 
(see below). Equation ([5]) implies that the modulus of detM is equal to 1. 

The entries of the transfer matrix are related to the reflection and transmission amplitudes of 
the potential V(x), denoted by r and t, respectively. When there are no waves coming from the 
right (B2 = 0), there exist solutions of the Schrodinger equation of the form 

e ik(x-a) + re -i*( a! -tt) j for x < fl) 

ip(x) = { te ik( - x ~ b \ for x > b, (6) 

u\ (x), for a < x < b, 

that is, solutions of the form ([2]) with A\ = 1, Ai = r, and B\ = t. Thus, from Eq. ([3]) it follows 
that 

1/i M 12 
r/t M 22 

with M\2 and M22 not yet identified, and from Eq. (JH) we obtain the well-known relation 



M 



1 



(7) 



(In most textbooks the reflection and transmission amplitudes are defined by means of expressions 
similar to Eq. ([6]), with r and t being the coefficients of e~ lkx and e lkx and therefore, the amplitudes 
r and t defined by Eq. (J6j) differ from those usually employed by factors e lka and e~ lkb , respectively.) 
Since V(x) is real, for k real the complex conjugate of the solution © 

r * e ik(x-a) _j_ e -ifc(a:-a)^ £ or x < a ^ 

t*e~ ik ( x - b \ for x > b, (8) 

u\(x), for a < x < b, 

is also a solution of the Schrodinger equation. Substituting the coefficients appearing in Eq. (JSj) 
into Eq. ([3]) we find that 

/ lit r*/t* 

M = 

{ r/t 1/t* 

Then, according to Eq. ([7]), det M = 1, which means that (for k real) the transfer matrix belongs 
to the group SU(1,1), formed by the 2x2 complex matrices with unit determinant that satisfy Eq. 

The reflection and transmission amplitudes of the potential V{x) for waves incident from the 
right, r' and t' , respectively, need not coincide with r and t. In fact, from Eqs. ([3|) and ([9]), setting 



(9) 



A\ = and B2 = 1, we must have 



\ 



I / l/t r*/t* 
I' I \ r/t l/t* 
which, making use of Eq. ([7|), implies that 

t' = t, ° = J + ¥ - ( 10 ) 
Hence, r = r' if and only if r/t is pure imaginary. 

3. Rectangular barriers 

The reflection and transmission amplitudes for a given potential are usually obtained by solving 
the time-independent Schrodinger equation (P) (see, e.g., Refs. 1-4). In the exceptional case of a 
piecewise constant potential, the transfer matrix (and, therefore, the reflection and transmission 
amplitudes) can be readily obtained by means of matrix exponentiation. 

The time-independent Schrodinger equation ([1]) can be expressed as the first-order differential 
equation 



_d_ 

d.r 




(11) 



where v(x) = 2mV(x)/fi 2 . Hence, if V(x) is a constant Vq for a < x < b, the solution of Eq. (lllj) is 
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i/j(x) 




" / 






j =exp 
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f{x) 









for a < x <b, where vq = 2mVo/h 2 , and c\, C2 are some constants. Thus, 

^ rur > 1 

cxp 



V 

with L = b — a. Letting 



V(a) 
V/(a) 



J 




(12) 



one finds that J 2 = (vq — fc 2 )J, where / is the unit 2x2 matrix, hence (see, e.g., Ref. 7) 



COsh(Ly/^]F) I - ^HLVW^) Jf i{vQ _ k 2 >Qi 

Vv - k 2 



exp(— LJ) = < 



it m \r sin(Wk' 2 - vq) t , 2 ^ _ 

cos(Lv — vq) I 7= , t — J, if f — k < 0, 



v 7 ^ 



(13) 



I -LJ, 



if v -k 2 = 0. 



On the other hand, from Eq. ([2]) we have 
^(a)=A 1 +A 2 , i;'(a) = ik(A 1 -A 2 ), 



^(b) =B X + B 2 , i>'{b) = \k(B x - B 2 ), 



that is, 





i -i/k Ma) m 

1 i/k ) { if/ {a) ) ' [ f{b) 

(Note that Eqs. (|14p correspond to the continuity conditions for tp and ij)' at x = a and x = b. 
Then, noting that 



(14) 



1 



-i/k 



J 



1 i/k 



1 1 

ik —ik 



1 
2ik 



v - 2k 2 

-Vq 



-V + 2k 2 



from Eqs. (|12p -(|14 p one finds that, for a rectangular potential barrier (or potential well) 



V{x) 



0, if x < a or x > b, 
Vq, if a < x < b, 



(15) 



the transfer matrix is given by 



M 



cosh(LVwo — k 2 ) I 



sin 



h(Wv - k 2 ) 1 



cos(L\/k 2 — vo) I 
/ , 



V^o - k 2 2ik 
sin(Lv / fc 2 — vq) 1 



VP 



2ik 




if v - k 2 > 0, 



if Vq - k 2 < 0, 



I 



\kh 



if Vq — k A 



0. 



(16) 



Note that owing to the definitions of the amplitudes A\, A 2 , Bi, and B 2 in terms of the exponentials 
e ±ik(x-a) anc | e ±ifc(a;— 6)^ transfer matrices (|16|) depend on a and b only through their difference 
L = b — a. The simplicity of the transfer matrices (|16p contrasts with the complexity of the 
expressions for the reflection and transmission amplitudes obtained in the standard manner (see, 
e.g., Ref. 2, chap. 5). Note also that even though we follow the conventions of Ref. 5, the transfer 
matrix (|16p does not agree with the amplitudes given in Eq. (12) of Ref. 5. 

It may also be noticed that, by allowing V^o — k 2 to become pure imaginary or taking the 
limit as \/vq — k 2 goes to zero, from the first expression in (|16p one can obtain the other two. 
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Furthermore, one can verify directly that, when k is real, the transfer matrices (|16p are of the form 
/ 

Cl.il 

, with \a\ 2 — |/3| 2 = 1 and therefore they indeed belong to SU(1,1). On the other hand, 



a (3 
(3* a* 

Eqs. (jlip ~ (|13p hold for k real or complex and since the trace of J is equal to zero for any value of 
k (even if V{x) was not real), the determinant of exp(— LJ) is equal to 1; therefore, the transfer 
matrices (|16|) have determinant equal to 1 also when k is complex, though M no longer belongs to 
SU(1,1). 

The example considered in this section also allows us to illustrate the fact that making use of 
the transfer matrix one can find the energies of the bound states or the quasinormal modes, by 
considering pure imaginary or complex values of k, respectively. 

In the case of the bound states of the potential well (|15p with Vo < 0, we have E < and 
writing k = i\k\, from Eq. ([2]) we see that in order for the wave function to remain bounded, B 2 = 
and A\ = 0. Then Eq. ([3]) implies that Mn, the first entry of the diagonal of M, must be equal 
to zero. Since the determinant of the transfer matrix is equal to 1 (independent of the value of k), 
this last condition is equivalent to saying that the off-diagonal entries of M (which have opposite 
signs) must be equal to +1 or —1. In the present case vq — k 2 < 0, and from the second line of Eq. 
([16]) we have 



sin(Lv/- k 2 - v ) v _ ±1 



tp(x) 



V-\k\' 2 ~vo 2\h 

which is equivalent to the conditions obtained in the textbooks (see, e.g., Refs. 1,2). 

The so-called quasinormal modes correspond to complex values of k for which there are no 
incident waves on the potential barrier but only outgoing waves. In this case the solution of the 
time-independent Schrodinger equation are of the form 

A 2 e- ik ^- a \ x<a, 

Bie ik(x-b)^ x>b ^ (17) 
u(x), a < x < b, 

assuming that the real part of k is positive (c/. Eq. ([2]), A\ and B2 are equal to zero so that there 

are no ingoing waves on the barrier). Thus, as in the case of the bound states, we have Mn = 

and, making use of the first expression in Eq. (|16p . we have 

v^r / H\ sinh^v^o - k' 2 ) v - 2k 2 

cosh(Lv^P) - v ^ _ k2 — = 

which can also be expressed in the form 



with the definition a = \/vq — k' 2 . Hence, a 2 + k 2 = vq. Following Chandrasekhar [8], we parame- 
terize k and a according to 

k = Qsina, a = Qcosa, (19) 

with Q 2 = vo and Q > (assuming vq > 0). Substituting these expressions for k and a into Eq. 
([TBI) we have 



. N Q 2 cos 2 a - Q 2 sin 2 a 

cosn(Lcr) — smh(-Lcr) —7, = U, 

2iQ z sin a cos a 

which is equivalent to 

cosh(Lcr) + i sinh(Lcr) cot 2a = 0. (20) 

Making use of the identities sinz = — isinh(iz), cos z = cosh(iz), this last equation can be written 
as sinh(Lcr) cosh(i2a) — cosh(Lcr) sinh(i2a) = 0, which is equivalent to 

sinh(Lcr — i2a) = 

and, therefore, La — \2a = in-zr, where n is an integer. Then, letting a = ot\ + i«2 and a = a\ + io"2, 
we have 

La\ = —2ct2, La% = 2a\ — mr. (21) 
From the relation a\ + \ai = Qcos(ai + \a<i) [see Eq. (|19p ] we obtain 

o\ = Q cos a\ cosh a%, 01 = — Q sin a± sinh a% (22) 

and, combining Eqs. (|2~T1) and (f22l) . it follows that 

— 2«2 = -^Q cos «i cosh «2> 2«i — n7r = —LQ sin ai sinh a<i. (23) 

Hence, 

2ai — nir 

tan ai tanh 02 = • (24) 

2a 2 

Similarly, from Eq. (fl~9l) . we have fci + i&2 = Qsin(ai + ia^), that is, 

k\ = Q sin a± cosh 02, k 2 = Q cos ot\ sinh 02 

and, making use of Eqs. (f23l) . 

7 2a2tanai 2a2tanha2 

fel = ^ , k 2 = - . (25) 
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By hypothesis, k\ > and Q > 0, therefore from Eqs. (|2"lj) and ([25]) it follows that 
if d2 > =>■ tanai < 0, 2a\ — nir < =>■ n > 0, — < ai < 7r, 

7T 

if a 2 < =>■ tanai > 0, 2a\ - nir > => n < 0, < a\ < — . 

Given a solution, ati, 02, of Eqs. (|23p . the values of fci and /C2 are determined by means of Eqs. 
(E5|). 

Since the time-independent Schrodinger equation ([T]) is obtained assuming that the wave func- 
tion has a time dependence of the form exp(— \Et/%). When k is complex, E has a negative 
imaginary part for k\ > [A;2 is negative, see Eq. (|25p ] that produces an exponential decay in time. 

Denoting by M^ a ' b ^ the matrix appearing in Eq. ([3]), we have the relation 

M (a,c) = M (a,b) M (b,c) ^ 

for any value of c. This relation together with Eq. (|16l) allow us to readily find the transfer matrix 
(or, equivalently, the transmission and reflection amplitudes) for any piecewise constant potential 
and from the condition M\\ = 0, the bound states and quasinormal modes can then be obtained, 
though the expressions will be even more involved than the ones considered here. 



4. Reflection and transmission of electromagnetic waves 



The behavior of a linearly polarized electromagnetic plane wave normally incident on a slab of 
dielectric material can be found following a procedure similar to that employed in the preceding 
section. For plane monochromatic waves propagating along the x-axis with the electric field parallel 
to the y-axis in a homogeneous dielectric medium, the wave equation reduces to 

d 2 K, 



dx 2 



+ k z E y = 0, 



(26) 



where k = nto/c, n is the refractive index of the medium and cu is the frequency of the wave. 
Equation ([26]) can be expressed as the first-order equation 

/ / , 1 \ / ; , ( 

(27) 



d 



En 



dx \ dEy/dx 



V 



1 

-k 2 



, dEy/dx 



which is of the form (jlip with v = 0; hence, 

/ 



E y (a) 
dE y /dx\ x=a 



M 



Ey(b) 

dE y /dx\ x=b 



where [see Eqs. CE} and dHJ] 



M = cos(fcL) I 



s'm(kL) 
k 



V 



1 

-k 2 



and L = b — a. 

If the slab is bounded by the planes x = a and x = b and, for instance, surrounded by vacuum, 
Eq. (|26p has solutions of the form 



Aie ik (x-a) + A 2 e- ik °( x - a \ for x < a, 

Bie ik Q (x-b) + g^-ikoix-b) for x > b 



(28) 



where k®= oo/c. Faraday's law imply that the z-component of the magnetic field is proportional to 
dE y /dx and, therefore, the continuity of the tangential components of the fields at the boundary 
of the slab amounts to the continuity of E y and dE y /dx and from Eq. ([28]) we see that 



E y (a) = A ± +A 2 , 



dEy 

dx 



(a) = ik (Ai-A 2 ), E y (b) = B 1 +B 2 , 



dEy 

dx 



(b) = ik (B 1 -B 2 ), 



thus, proceeding as in the previous section we obtain the relation 




M 



Si 



B 2 



(29) 



M = cos(fcL) I 



s'm(kL) i 
k Zk Q 



( 



with the transfer matrix 

k 2 + k 2 -k 2 +k 2 \ 

U2 _ 1.2 u2 _ r.2 I 

which is related to the transmission and reflection amplitudes as in Eq. ©. It may be noticed 



(30) 



that, also in the present case, the transfer matrix (|30p belongs to SU(1,1) for k real and that 
the determinant of M is equal to 1 even if k is complex (which would correspond to a nonzero 
conductivity) . 



Acknowledgment 

One of the authors (I.R.G.) thanks the Vicerrectorfa de Investigacion y Estudios de Posgrado of 
the Universidad Autonoma de Puebla for financial support through the programme "La ciencia en 
tus manos." 



10 



References 



1. L.I. Schiff, Quantum Mechanics, 3rd ed. (McGraw-Hill, New York, 1968). 

2. S. Gasiorowicz, Quantum Physics, (Wiley, New York, 1974). 

3. I.I. Gol'dman and V.D. Krivchenkov, Problems in Quantum Mechanics, (Pergamon, London, 
1961; reprinted by Dover, New York, 1993). 

4. K. Gottfried and T-M. Yan, Quantum Mechanics: Fundamentals, 2nd ed., (Springer, New 
York, 2003). 

5. L.L. Sanchez-Soto, J.F. Carinena, A.G. Barriuso and J.J. Monzon, Eur. J. Phys. 26 (2005) 
469. 

6. J.R. Reitz, F.J. Milford and R.W. Christy, Foundations of Electromagnetic Theory, 4th ed. 
(Addison- Wesley, Reading, MA, 1993). 

7. K.B. Wolf, Rev. Mex. Fis. 49 (2003) 465. 

8. S. Chandrasekhar, Proc. R. Soc. London A 344 (1975) 441, reprinted in Selected Papers, Vol. 
6, (Chicago University Press, Chicago, 1991). 



11 



